MIT-CTP/4127 



A controlled expansion for certain non-Fermi liquid metals 
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The destruction of Fermi liquid behavior when a gapless Fermi surface is coupled to a fluctuating 
gapless boson field is studied theoretically. This problem arises in a number of different contexts 
in quantum many body physics. Examples include fermions coupled to a fluctuating transverse 
gauge field pertinent to quantum spin liquid Mott insulators, and quantum critical metals near 
a Pomeranchuk transition. We develop a new controlled theoretical approach to determining the 
low energy physics. Our approach relies on combining an expansion in the inverse number (TV) of 
fermion species with a further expansion in the parameter e = zt, — 2 where Zb is the dynamical critical 
exponent of the boson field. We show how this limit allows a systematic calculation of the universal 
low energy physics of these problems. The method is illustrated by studying spinon fermi surface 
spin liquids, and a quantum critical metal at a second order electronic nematic phase transition. 
We calculate the low energy single particle spectra, and various interesting two particle correlation 
functions. In some cases deviations from the popular Random Phase Approximation results are 
found. Some of the same universal singularities are also calculated to leading non- vanishing order 
using a perturbative renormalization group calculation at small TV extending previous results of 
Nayak and Wilczek. Implications for quantum spin liquids, and for Pomeranchuk transitions are 
discussed. For quantum critical metals at a nematic transition we show that the tunneling density 
of states has a power law suppression at low energies. 

PACS numbers: 



I. INTRODUCTION 

This paper is concerned with the destruction of Fermi 
liquid behavior in two dimensional systems where a gap- 
less Fermi surface is coupled to a fluctuating gapless bo- 
son field. This problem arises in a number of differ- 
ent contexts in quantum many body physics. A well 
known example is where the fermions are coupled to a 
gapless transverse gauge boson. This describes the low 
energy effective theory of certain quantum spin liquid 
phases^, the theory of the half-filled Landau level^, and 
various non-fermi liquid metallic phases^— . A different 
and equally well known example is as a description of 
quantum critical metals at a 'Pomeranchuk' instability. 
The classic example is the Stoner transition associated 
with the onset of ferromagnetism in a metal. In recent 
years attention has focused on a different example of a 
Pomeranchuk transition: that associated with the onset 
of electronic nematic order— ~— from a Fermi liquid metal. 
Here electronic nematic order means a phase where the 
lattice point group symmetry but not translation sym- 
metry is broken. Such order has been observed with in- 
creasing frequency in a number of different correlated 
metals^— giving rise to an interest in the associated 
quantum phase transition. At such a quantum phase 
transition the nematic order parameter is described as a 
gapless fluctuating Bose field, and its coupling to the gap- 
less Fermi surface destroys Fermi liquid behavio r 10 i 12 ' 13 . 

The purpose of this paper is to formulate a new con- 
trolled theoretical approach to this class of problem, the 
need for which has been emphasized recentl y—^" . The 
low energy physics of the resulting non-fermi liquid metal 
is characterized by universal scale invariant behavior. 
Our approach provides a systematic method of calculat- 



ing the exponents and other universal properties asso- 
ciated with this scale invariant behavior. We illustrate 
this by studying many physical properties of the gauge 
field model and of the nematic quantum critical metal in 
detail. 

Quite generally the low energy physics of problems of 
this sort is conveniently described by restricting attention 
to fermionic modes in the immediate vicinity of the Fermi 
surface, and the long wavelength, low frequency modes 
of the fluctuating boson field. The model is described by 
the Euclidean action 
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Here fk a , & — 1,....N, is a fermion field with N possible 
flavors and a is the boson field. In the gauge model a 
is the transverse component of a U(l) gauge field, and 
0(x, t) is the transverse component of the current den- 
sity of fermions. At a nematic quantum critical point a 
will be taken to be the nematic order parameter field, 
and 0{k,uj) is the fermion bilinear with the same sym- 
metry. For instance, on a two dimensional square lattice 
with lattice constant I a uniform nematic order parame- 
ter couples to N^ 1 / 2 J2 k (cos(k x £) - cos(k y £)) fka fka- 

Much prior work of course exists on this problem. In a 
number of early papers&^ir— the problem was analysed 
in a Random Phase Approximation (RPA) and various 
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related approaches. This showed that the fermions and 
gauge bosons stay strongly coupled in the low energy 
limit. In the RPA, the boson propagator is overdamped 
due to Landau damping by the gapless Fermi surface. 
The fermion self energy has a power law non-fermi liq- 
uid frequency dependence. Further the long wavelength 
density response function retains its Fermi liquid form. 

In the gauge field problem, some of the RPA results 
were further substantiated 27 through a quantum Boltz- 
mann approach which considered the fate of various pos- 
sible shape fluctuations of the Fermi surface. Smooth 
shape deformations of the Fermi surface (which deter- 
mine long wavelength density response and the gauge 
propagator) were shown to retain Fermi liquid behav- 
ior while "rough" deformations have the potential to be 
non-fermi liquid like. The latter determine the behavior 
of the single fermion Green's function and the structure 
of the 2Kf singularities (i.e at wavevectors connecting 
antipodal tangential portions of the Fermi surface) in re- 
sponse functions. These main results were further sup- 
ported in comprehensive diagrammatic analyses 28 of the 
model which suggested that the leading RPA answers 
for many quantities were in fact exact in the low energy 
limit. In particular the structure of the gauge field prop- 
agator, the fermion self energy, and the long wavelength 
density response were argued to have the same form as 
the RPA result. Similar diagrammatic analyses with the 
same conclusions have also been reached 13 for the ne- 
matic quantum critical point. In the gauge held case the 
2Kf singularities in the density response function were 
argued to have specific non-fermi liquid like power law 
forms 28 . 

Is there a controlled limit in which the reliability 
of these results may be assessed? One attempted 
approac h 28 ! 29 is to take the limit of N (the number of 
fermion species) large, and expand in powers of 1/N. In 
the early work it was argued 29 that at low energies in 
the large- N limit only patches of the Fermi surface with 
parallel normals are strongly coupled to each other. Any 
such patch couples strongly to a boson whose momen- 
tum is perpendicular to the normal to the Fermi surface. 
The low energy physics is therefore correctly described 
by focusing attention on patches with parallel normals. 

In some remarkable recent work, Sung-Sik Lee^ 2 - reex- 
amined the model of N fermion species coupled to a U(l) 
gauge field in the large- N limit. He showed that even at 
large- TV the theory remains strongly coupled, and that 
its solution requires non-trivial summation of an infinite 
number of Feynman diagrams. When only a single patch 
of the Fermi surface is considered, a book-keeping device 
was introduced to show that the 1/N expansion could be 
organized in terms of the genus of the surface in which 
the Feynman diagrams were drawn. Based on this the 
general validity of the physical picture built up by RPA 
and the other earlier analyses for small N has been ques- 
tioned. 

Even more recently in another very interesting paper 
Metlitski and Sachdev 23 studied the fate of the theory 



with both a Fermi surface patch and its antipodal part- 
ner included. This is believed to be fully sufficient to 
correctly describe the asymptotic low energy physics of 
the system. They found a number of further difficulties 
with the large- N expansion. Specifically higher loop cor- 
rections for the gauge propagator involved higher powers 
of N than the leading order one loop RPA result. This 
unpleasant finding led them to question the existence of 
a well-defined large- TV limit to control the theory. These 
authors also showed that in a perturbative loop expan- 
sion the self energy acquires singular momentum depen- 
dence at three loop order. However the loop expansion 
has no apparent control parameter. 

In light of these results it becomes important to search 
for alternate reliable methods to judge the validity of 
RPA and other diagrammatic approaches to the problem. 
In this paper we introduce a new controlled expansion 
to determine the low energy physics of this model. We 
consider a family of models where the 'bare' boson action 
is modified to 

r |?| z f>- 1 

L L± ^Hk,u J )\ 2 (5) 

J k,u e 

The number Zf, (the boson "dynamical critical expo- 
nent") equals 3 in the original model in Eqn.Q] The case 
Zb = 2 arises in the theory of the half-filled Landau level 
with long range 1/r Coulomb interactions between the 
electrons^, and in the theory of the bandwidth controlled 
Mott transition of the half-filled Hubbard model devel- 
oped in Ref. H3- We show that the large- N expansion 
can be controlled in the limit of small e = 2% — 2. Specif- 
ically we show that the limit N — > oo,e = Zf, — 2^0 
such that eN is finite leads to reliable answers for the 
low energy behavior of the system. We demonstrate that 
the RPA answers for the fermion and boson propagators 
are indeed exact in this limit. A systematic expansion in 
powers of 1/N is possible for small e. Deviations from 
RPA emerge at higher orders in the 1 /N expansion. Fur- 
thermore differences between the gauge model and the 
nematic critical point also appear. At order 1/N 2 , we 
find a singular momentum dependent correction to the 
fermion self energy - however, in the gauge model, this 
singularity is subdominant to the leading order momen- 
tum dependence, so that the fermion self energy retains 
its RPA form, at least to this order. Further we calculate 
the exponent characterizing 2Kf particle-hole singulari- 
ties, and show in the gauge model that (depending on 
the value of eN), they may be enhanced compared to a 
Fermi liquid. For a quantum critical metal at a nematic 
transition the fermion propagator again retains its RPA 
form at leading order but at o(l/N 2 ) acquires a singu- 
lar correction to the self energy that dominates over the 
RPA form. This modification from the RPA is in ac- 
cord with the calculation of Ref. [23| but is now performed 
in a controlled expansion. A further difference with the 
gauge field problem is in the structure of the 2K f singu- 
larities. We present calculations and physical arguments 
that show that the 2Kf singularities are weakened at the 
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nematic quantum critical point compared with a Fermi 
liquid. 

A crucial physical ingredient that determines the low 
energy physics is the nature of the 'Amperean' interac- 
tion between two fermions that is mediated by the boson 
field. The term 'amperean' is appropriate for the gauge 
field case where the interaction is between fermion cur- 
rents but we will use it to describe the nematic transition 
as well. In the gauge field case, the currents of a particle 
in one patch of the Fermi surface is parallel to that of a 
hole in the antipodal patch. By Ampere's law the gauge 
mediated interaction between such a u 2Kf" particle-hole 
pair is attractive. In contrast the currents of two particles 
with one from either patch are antiparallel and the gauge 
mediated interaction is repulsive in the particle-particle 
(Cooper) channel. The situation is reversed in the ne- 
matic model. We show how this difference between the 
particle-hole and Cooper channel interactions plays an 
important role in many aspects of the low energy physics. 

Previously Nayak and Wilczek 31 - studied the low en- 
ergy physics of the gauge field model for small e, and 
finite N using perturbative renormalization group (RG) 
methods. The one loop beta function for the coupling 
constant e takes the form 

/3(e 2 ) = f e 2 - ^e 4 (6) 

where c is a positive constant. This leads to a perturba- 
tively accessible non-fermi liquid fixed point for iVe small 
and positive. Some properties of this fixed point were cal- 
culated in Ref. l3l| and shown to be consistent with the 
RPA analysis. Calculations with this one loop beta func- 
tion are able to provide answers for the scaling exponents 
to order e for any N. In contrast our approach of directly 
solving the theory at large- TV enables us to extract expo- 
nents that are high order or even non-perturbative in e. 
Thus for instance the deviations from RPA (subdominant 
in the gauge field model) discussed above are expected 
to appear in the e expansion only at order e 3 , and so do 
not show up in the o(e) calculations. In the regime where 
they overlap (i.e to order e at fixed large- N), we will show 
that the exponent values for many properties calculated 
within our approach agree with those obtained from the 
perturbative RG. As it does not seem to be available 
in the literature, we calculate the 2Kf exponent at this 
small e, finite N fixed point. We show that the exponent 
has an interesting non-analytic dependence on e for small 
e which leads in the gauge model to an enhancement of 
the 2Kf singularities compared with a Fermi liquid in 
this limit. In the nematic model we find a suppression 
of the 2Kf correlations as expected on general physical 
grounds that we also discuss. 

An important feature of the low energy physics is that 
the Fermi surface is preserved and is sharp even though 
the Landau quasiparticle is destroyed. Further at low 
energies and for momenta close to the Fermi surface the 
fermionic spectrum is scale invariant. This was already 
implied by the RPA results, and survives in our treate- 
ment. A similar picture was also argued 3 ^ to describe 



continuous phase transitions where an entire Fermi sur- 
face disappears (such as a continuous Mott transition). 
Following Ref [32| we will refer to this as a critical Fermi 
surface. Some (though not all) aspects of critical fermi 
surfaces associated with Mott-like transitions may be ex- 
pected to be shared with the Pomeranchuk transitions 
discussed in this paper. It is therefore useful to consider 
these results in terms of a general scaling form expected 
for fermions with a 'critical Fermi surface'. We write for 
the fermion Green's function 

G <*->~Rf'°(f) < 7 > 

Here kn is the deviation of the momentum from the Fermi 
surface. Note that the z that enters this scaling equation 
is the 'fermionic' dynamical critical exponent. For the 
problems studied in this paper the RPA gives z — 4j and 
a = 1 . The latter is a result of the absence of any singular 
momentum dependence in the self energy in RPA. Our 
results may be viewed as a calculation of a and z within 
a systematic expansion. In the gauge model to o(l/N 2 ) 
these exponents do not change as far as the leading singu- 
lar structure is concerned. For the nematic critical point 
we find that a = 1 — r\f with r/f positive. 

Following the general discussion in Ref. HH, we show 
that the difference from the RPA result has direct and 
measurable consequences for the electron single particle 
tunneling density of states N(E) (where E is measured 
from the chemical potential) at the nematic quantum 
critical point. Within the RPA the tunneling density of 
states is a constant at the Fermi level. However beyond 
RPA there is a power law suppression of N(oj): 

N(w)~\wft (8) 

The exponent rjf is calculated in Section [IVI Extrapola- 
tion of the leading order results to z\, = 3, N = 2 gives 
the estimate r\f f=a 0.3. 

What about the fate of the large- N limit when Zb is 
not close to two? We suggest that recent calculations of 
Ref. should be interpreted as an instability towards 
translational (and possibly other) symmetry breaking in 
this limit. In the nematic context this means that there 
is no direct second order quantum phase transition as- 
sociated with nematic ordering in two dimensions if N 
is sufficiently large. Instead the transition is preempted 
by the appearance of density wave and possibly other 
orders. So if for the physical case N — 2 there is a di- 
rect nematic transition, then it is not usefully accessed 
by the large- N expansion. Our approach of combining 
the large- N with a small Zb — 2 or a direct small Zb — 2 
perturbative RG then become the only available meth- 
ods to theoretically access such a quantum critical point 
directly in two dimensions. Similar phenomena also hap- 
pen in the gauge field model - in the large-iV limit we 
propose that the uniform state is unstable to translation 
symmetry breaking. 
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FIG. 1: (color online) Our suggested phase diagram. Above 
the indicated curve, the putative critical theory is likely pre- 
empted by some other broken-symmetry state. The behavior 
of the proposed phase boundary at small N is one possible 
extrapolation. 




FIG. 2: (color online) An alternate possible phase diagram. 
Here the interesting point zj, = 3, N = 2 is in the unstable 
regime; in this case, it would be best accessed starting from 
the correct mean field theory for the new broken-symmetry 
state. 



It is instructive to consider the behavior of the model 
in a two dimensional plane spanned by Zb and 1/N. We 
show our proposed 'phase diagram' in Figs. Q] and [21 It 
is clear that the approach developed in this paper is ide- 
ally suited to describing the gauge model or the nematic 
quantum critical point for N = 2 if it is not part of the 
unstable region, i. e if Fig. Q] is realized. If on the other 
hand Zb — 3, N = 2 belongs to the unstable region as 
depicted in Fig. [21 we may still hope that the approach 
in this paper is useful in describing the physics at tem- 
peratures above the onset of the instability. 

The rest of the paper is organized as follows. In Sec- 
tion. [TTl we begin with some preliminaries and briefly 
discuss the patch construction for the Fermi surface that 
is used in the rest of the paper. Some subtle but im- 
portant aspects of the patch construction are relegated 
to Appendix [XJ Then in Section Mil we warm up by 
studying the theory of just one patch and ignoring any 
coupling with the other antipodal patch, and show how 
our expansion provides a controlled answer in this sim- 
plified problem. We then study the full two patch the- 



ory in Section IIVI and determine the singular structure 
of the boson and fermion propagators. In Section [V] we 
present a calculation of the exponent characterizing 2Kf 
singularities within our approach. Next in Section [VTl we 
explore the connections with the perturbative RG calcu- 
lations of Ref. and extend their results to 2Kf singu- 
larities. In Section IVIII we discuss simple physical inter- 
pretation of the results of the calculations and their con- 
sequences. Section IVIIII describes our suggestions on a 
possible phase diagram. We conclude in Section HXl with 
a general discussion on how our results fit in with var- 
ious other related problems and theoretical descriptions 
of non-fermi liquid metals. Various appendices contain 
details of calculations. 



II. PRELIMINARIES 

As mentioned above the low energy physics is cor- 
rectly described by focusing attention on Fermi surface 
patches with parallel normal o 8 ' 23 i 28 ' 29 . This is because 
the interactions mediated by the boson field are predom- 
inantly small-angle scattering processes. Furthermore 
short range four fermion interactions that couple differ- 
ent patches become unimportant at low energie o 8 ' 28 ^ 29 as 
can be checked a posteriori after the two patch theory is 
solved. Thus the universal low energy physics of the sys- 
tem is correctly captured by a theory that focuses on two 
opposite patches of the Fermi surface. We focus hence- 
forth on two opposite patches of Fermi surface; there are 
a number of subtle and important points related to the 
patch construction that we elaborate on in Appendix [A] 
The patch construction also has a number of immediate 
consequences for the behavior of many physical proper- 
ties. These will be discussed in Section IVlTl 

Consider patches of the Fermi surface with normals 
along ±x. We will denote the corresponding fermion 
fields fn/L where R denotes the right patch and L the 
left one. It is useful to begin by considering the boson 
and fermion propagators in perturbation theory keeping 
just the leading one loop diagrams (Figs. [3] and The 
imaginary frequency boson propagator D(k,u) becomes 

D(k,uj) = — — — — (9) 

T|fc„| + 

with 1 7 = t-. Unless otherwise mentioned we will hence- 
forth set e = 1. The fermion propagator is determined 



For the case of the nematic, we absorb the dependence on the 
angle between the nematic ordering vector and the patch in ques- 
tion into the coupling e. This coupling e specifies an energy scale 
below which our low energy theory is applicable. Since this en- 
ergy scale vanishes at the 'cold spots', we must restrict attention 
to a patch of the Fermi surface away from this direction. 
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by its self energy which at one-loop level takes the form 

E = -t-^sgn(w)H^ (10) 
The constant A is given by 

A = 47rsin — 7 »t (11) 

Zb 

and thus vanishes linearly as zi, — > 2. In terms of the scal- 
ing form in Eqn. [7] this implies the fermionic dynamical 
critical exponent z — 4^ and a — 1 as promised. 
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FIG. 3: 1-loop boson self-energy. 
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FIG. 4: 1-loop fermion self-energy. 

The arguments of Ref. [H show that a minimal Eu- 
clidean action that enables correct description of the low 
energy physics is given by 

S = S f + S mt + S a (12) 
S f = d 2 xdr ^2 fsa {rid T - isd x - d 2 ) f sa (13) 

J sa 

Sint = / d 2 xdT—^=af sa f sa (14) 
J v N 

Sa = I Iky^-^afcu,)] 2 (15) 

Here s — +1 for the patch R and —1 for the patch L. The 
parameter rj is taken to be small and positive. The field a 
represents just the x-component of the vector field at. In- 
deed it is just this component that couples strongly to the 
patches with normals along ±x. Note in particular that 
the boson field couples with opposite sign to the two an- 
tipodal patches. If on the other hand we were interested 
in the critical theory for a Pomeranchuk transition (such 
as a transition to an electronic d-wave nematic state in a 
two dimensional metal which microscopically has square 
lattice symmetry), the minimal action will have a very 
similar form except that the boson will couple with the 
same sign to antipodal patches. While this difference is 
unimportant for some properties it plays a crucial role 
in others. For instance the structure of the 2Kf singu- 
larities is completely altered between the gauge field and 
nematic models. 



III. ONE PATCH THEORY 

We begin by focusing attention only on one patch, say 
the right one, and completely ignoring the other one. In- 
deed Ref. [12 showed that the standard large- A expansion 
leads to an apparently strongly coupled theory already in 
this simplified model. The main point is that a high loop 
diagram may formally look like it is high order in the 
1/A expansion. However for many such diagrams the 
corresponding loop integral diverges in the rj — > limit. 
This divergence may be regularized by using the one loop 
self energy in the fermion propagator. As this is of or- 
der 1/N, the singular r\ dependence is traded for an en- 
hanced power of A in the numerator. Consequently the 
naive 1/A counting is modified and an infinite number 
of diagrams survive in each order of 1/N. A systematic 
way to keep track of the true power of 1/N is obtained by 
using a "double-line" representation for the boson field 
that was previously used in the treatment of the electron- 
phonon interaction in metal s 33 ' 34 . It was shown that the 
1/N expansion could be organized as a genus expansion 
with all "planar" diagrams surviving to leading order. 
Ref. further established that in the large- N limit the 
boson propagator is unrenormalized beyond 1-loop - in 
other words all higher loop diagrams that survive in the 
large-A limit give vanishing contributions. Each indi- 
vidual term contributing to the fermion self energy is (if 
one calculates using the 1-loop fermion propagator) fi- 
nite, and has the same functional form as the 1-loop self 
energy: formally (at Zb = 3), 

S = -il S gn(a;)Ml^^^ (16) 

n 

The nth term in the sum comes from diagrams that are 
formally of order 1/N n in the large- N expansion. How- 
ever for all planar diagrams there is a compensating en- 
hancement factor A" -1 in the numerator so that each 
term is of order 1/N. The worry is whether the sum over 
the infinite contributing diagrams leads to something sin- 
gular or not. 

It is straightforward to see that these results carry over 
to general Zb- Indeed the kinematics leading to the diver- 
gences in the small n limit depend only on the existence 
of the gapless Fermi surface and not on the detailed form 
of the boson propagator. When the divergence is regu- 
larized with the one loop fermion self energy, every — is 
traded for a factor AA. Eqn. [11] is accordingly modified 
to 

1 . ,. 2 b n (AA)" -1 
E = -i-sgn(w)M'»£-^^ (17) 

n 

Here the coefficients b n are all independent of A but in 
general depend on 2%. The utility of the small Zb — 2 
limit where A oc Zb — 2 is now apparent. So long as Zb — 2 
is of order 1/A, the enhancement factor (AA) ra_1 in the 
numerator of each term above is finite. If further the 
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6„'s are sufficiently non-singular when Zb — >• 2 then in 
the large- iV limit only the n = 1 term survives. This is 
just the one loop answer which is thus exact in this limit. 
This claim can be illustrated explicitly by calculating a 
particular instance of a dangerous diagram which con- 
tributes to the series above at n = 2, such as the one 
shown in Fig. [T3] We do this in Appendix [B] and show 
that though it is of order l/N for general Zb, in the limit 
Zb — 2 cx l/N it becomes higher order in l/N. In partic- 
ular the corresponding coefficient 62 has a finite limit as 
Zf, — > 2 so that its contribution is of order l/N 2 . This is 
in fact expected to be true for all the b n which have limits 
when Zb — > 2 such that to leading order high n terms in 
the series give subdominant powers of N to the leading 
order result. Indeed the absence of wlog 2 u; terms ^ 28 i 35 
in the self energy exactly at Zb = 2 implies that the self 
energy for small A at fixed N has at most one inverse 
power of A. Furthermore in this limit the self energy can 
be explicitly calculated using a perturbative RG tech- 
nique (see Section IVTf . The answer agrees exactly with 
the leading order term in the series above. This means 
that the small Zb — 2 behavior of b n is such as to keep 
the high n terms subdominant to the leading order one 
in the large- TV limit. 

The boson propagator is also given exactly by the one- 
loop answer. For arbitrary Zb this follows from the ar- 
guments of Ref. [22] due to the vanishing of higher order 
planar diagrams. However in the small Zb — 2 limit it also 
follows for the same reason as above - the enhancement 
factors that render higher order diagrams to be of the 
same nominal order in 1 /N all become innocuous in the 
small Zb — 2 limit. 

We emphasize that even though we have used the for- 
mal device of small Zb — 2 to control the large- N expan- 
sion the frequency dependence of the fermion self energy 
2 

£ oc — isgn(u;)|u;| z b is exact to all orders in z& — 2. The 
smallness of Zb — 2 merely assures us that the proportion- 
ality constant has a sensible l/N expansion. 

Similar results also apply to the fermion-boson vertex 
which to leading order in 1/7V is unrenormalized. 

It is useful to understand the scaling structure of the 
low energy physics described above. We have shown that 
in the one patch theory the large- N, small Zb — 2 limit, the 
low energy physics is described by a fixed point invariant 
under the scaling transformation 
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(20) 
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(21) 


J) 




(22) 



This is the same scaling structure that is obtained in a 
naive one loop approximation. Though Ref. [22| has raised 
concerns over whether this scaling is internally consistent 
in the large- N limit we see from the preceding analysis 



that it indeed is if the limit of small Zb — 2 is also simul- 
taneously taken. 

The considerations of Ref. [28| can now be used to argue 
that this scaling structure is exact to all orders in the 
1 /N expansion (so long as Zb — 2 is small) within this one 
patch theory. More specifically the boson propagator and 
fermion self energy will have the same functional forms 
(Eqn. [5] and Eqn. [TU)) as the one loop answers. Thus we 
expect that in the N, Zb plane there is a region of finite 
extent where the scaling structure above is preserved. 
This conclusion is further bolstered by the perturbative 
RG analysis for finite N, small Zb — 2 (see Section IVTi) . 



IV. TWO PATCHES 

The considerations above are readily generalized to the 
two patch theory. Indeed as before at each order of l/N 
there are a number of diagrams that are divergent in the 
small 77 limit. This is traded for an enhancement fac- 
tor of AiV raised to some power when the one-loop self 
energy is used instead. For Zb — 2 of order l/N these 
enhancement factors become finite and a controlled l/N 
expansion emerges. The structure of the leading order 
contribution in 1 /N to the fermion self energy, the boson 
propagator, or the interaction vertex is then not modi- 
fied from the one patch theory, and thus retains its RPA 
form. In particular the scaling structure in Eqns. IT8"ll2"2l 
is preserved. Below we examine higher order corrections 
in the l/N expansion, and show that new singularities 
appear at o(l/N 2 ). These could potentially modify the 
scaling structure from that in Eqns. [TSll2"2l above. How- 
ever we show that in the gauge field model these con- 
tribute only to subdominant corrections to the one loop 
fermion Green's function. On the other hand for the 
nematic transition there is indeed a modification of the 
fermion scaling 'dimension' at o(l/N 2 ). Our calculations 
rely on the very recent impressive results of Metlitski and 
Sachdev 23 who studied the boson propagator and fermion 
self energy in a direct (albeit uncontrolled) perturbative 
loop expansion upto three loop order. Here we will show 
that calculations along the lines of those in Ref. [23J leads 
to controlled results for various physical quantities within 
our modified l/N expansion. 

First consider the boson propagator. At Zb — 3 the 
analysis of Ref. established that to three loop order 
there is no shift of the true dynamical critical exponent. 
This was done by showing that the boson propagator at 
zero external frequency stays proportional to q 2 up to 
three loops. However the three loop contribution is of 
order y/~N bigger than the one loop contribution thereby 
casting doubts on the existence of a sensible large- iV 
limit. We will suggest an interpretation of this result 
in Section IVIIII But for now we discuss the results of an 
identical analysis in our limit of small Zb — 2 ~ o(l/N). 
First we note that the leading order term in the inverse 
gauge propagator is of order 1 with our conventions. Next 
higher loop diagrams clearly give subdominant powers of 
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1/N as there are no enhancement factors in the limit of 
small Zfy — 2. At zero external frequency the leading 1/N 
correction comes from the two diagrams shown in Fig. [5] 
These same diagrams were calculated in Ref.HHat Zb = 3. 
Repeating for general Zf,, we find the 1/N correction 



±c 



A 1 



(23) 



with the — sign for the nematic critical point and the + 
sign for the gauge model. The function /i is evaluated 
in Appendix [C] , and is readily seen to have a finite limit 

when z b -> 2. For large XN, we have fi(XN) ~ (AA)^" 
in agreement with the result of Ref. [23| when Zb = 3. We 
see explicitly that when z\, — 2 is o(l/N) these three loop 
contributions are down by a factor 1 /A compared to the 
one loop term. Thus the large-A expansion is indeed well 
defined in this limit. 



! X 





FIG. 5: Three-loop boson self-energy diagrams. 



FIG. 6: Two-loop fermion self-energy diagrams merely renor- 
malize the coefficient of loil 2 ^ 26 . 



P + ^ £ 1 



k + q k + l 



k-q 



FIG. 7: Three-loop Fermion self-energy diagrams involving 
fermions on both patches. Fermions on the right patch are 
denoted by solid lines, fermions on the left patch by dashed 
lines. 



(for the right moving fermion) 



5V(p,u> = 0) - ± — J(\N)( Px +p 2 y )\n 



where the function J(AA) is defined in Appendix lO and 
is positive definite. The + sign applies to the nematic 
critical point and the — sign to the gauge model. In a 
renormalization group framework this can be interpreted 
as the leading terms of a singular contribution to the self 
energy of the form 




(25) 



A similar contribution exists in the frequency dependent 
part as well, consistent with the dynamical scaling. In 
the gauge field case the plus sign applies, and this singu- 
lar correction is subdominant to the terms that already 
exist at leading order. Indeed we expect that any per- 
turbation of the fixed point by irrelevant perturbations 
will generate an analytic contribution to the momentum 
dependence of the self energy that will then dominate 
over the singular corrections found at o(I/A 2 ). Though 
the leading order frequency dependence is not analytic 
we expect that if we use the large-A fermion propaga- 
tors to calculate the effects of an irrelevant operator in 

2 

perturbation theory we will simply again generate a |w| '>> 
term. This will dominate over the singular order I/A 2 
correction. Thus we conclude that in the gauge model 
the leading singularities are correctly given by the RPA 
forms at least to order 1/N 2 . We note that our interpre- 
tation is different from that in Ref. [23l 

In the nematic case the minus sign applies in the ex- 
ponent of Eqn. 1251 This is more singular than the 'bare' 
momentum dependence of the inverse Green's function, 
and consequently will dominate the low energy physics 
near the Fermi surface. This singular correction can be 
interpreted as a shift of the scaling of the fermion fields 
from that in Eqns. [T%ll2"2l Thus to order 1/A 2 we have 



2T fa +5-T7 f 

fLiPxiPy,"') = b 4 f a s(px,Py,u) 



(26) 



where 



Next we consider the fermion propagator. The one- 
loop self energy is inversely proportional to A A and hence 
is of order I. It is easy to see by explicit calculation that 
the two loop diagram shown in Fig. [6] is momentum inde- 
pendent, and merely provides an o(l/A 2 ) modification of 
the coefficient of the frequency dependent part of the self 
energy. The most important effect at this order comes 
from the two graphs shown in Fig. [7] At zero exter- 
nal frequency these graphs lead to singular momentum 
dependence 2 ^,. Details are given in Appendix [C] We find 



3A : 



rJ(AA) 



(27) 



All the other scaling equations remain unmodified. This 
implies that the fermion Green's function satisfies the 
scaling form with a = 1 — rj t . 

What is the physical origin of the signs and the dif- 
ferences between the gauge and nematic models? We 
explain this in Section IVII1 To set the stage we first cal- 
culate singularities in some other quantities within the 
general two patch theory. 
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V. 2Kf AND OTHER SINGULARITIES 

The calculations of Ref. [28| on the gauge model showed 
that the response to an external field that couples to the 
fermion density at momentum 2K t is modified from that 
of a Fermi liquid due to the gauge interaction. The phys- 
ical origin of this effect is clear. First there is a suppres- 
sion of the 2Kf response coming from the smearing of 
the Landau quasiparticle due to the boson interaction. 
Second there is (in the gauge model) an enhancement 
coming from the "Amperean" attraction between a par- 
ticle at Kf and a hole at —Kf. The gauge currents of 
such a particle and hole are parallel to each other so that 
the gauge interaction is attractive for such a particle- 
hole pair. The net modification of the 2Kf singularity 
is determined by the interplay between these two effects. 
We emphasize that the Amperean enhancement is spe- 
cific to the gauge field problem. For the closely analo- 
gous problem of a quantum critical point associated with 
a Pomcranchuk transition, the particle-hole interaction 
mediated by order parameter fluctuations is repulsive. 
This goes in the same direction as the effect due to the 
smearing of the Landau quasiparticle. So there is no 
competition and we expect that the 2Kf singularities are 
simply suppressed when compared with the Fermi liquid 
at the Pomeranchuk transition. 

The structure of the fixed point in the large- N, small 
Zb — 2 limit enables a controlled calculation of these ef- 
fects. Consider an external field that couples to the 2Kf 
fermion density through the following term in the action 



d xdrf L f R + h.c 



(28) 



By power counting it is easy to see that the scaling in 
Eqns. H8ll22l implies that the coupling u scales as 



i = ub 



(29) 



In general this will be modified at order l/N as we 
demonstrate below. But first let us understand how the 
scaling of the u determines the 2Kf singularities. Assume 
in general that u scales as 



This implies that the operator P2K f (x, t) 
j 'l(x ', t) j 'r(x ', t) scales as 



f>2K f ( x ' \v' V) = b A P2k f (x,y,T), 



with 



A 



Zb 



(30) 



(31) 



(32) 



This determines the behavior of the singular part of 
the 2Kf density correlation function CiKj (x, y, t) = 
(p*2K s { x > Vi T )P2K f (0, 0, 0)). Its Fourier transform satis- 
fies 



C'2K f (Px,Py,u) 



- 2A CU f {p' x ,p' y ,w') (33) 



Note that the momenta p x ,Py describe the deviation of 
the full momentum from 2KfX in this correlation func- 
tion. We then immediately have the scaling form 



C2K f {Px,Py,u) 



-C 



\P( 



r" pi 



Note that in the usual Fermi liquid case, we have Zb = 2, 
4> u = 1 which reproduces the well known square root 
frequency dependence of the singular part of the 2 Kf 
correlations. 

The leading order correction to <j) u in the l/N expan- 
sion comes from the one loop vertex correction diagram 
of Fig. [15] To calculate it we combine the large- TV expan- 
sion with an RG transformation in which internal loop 
integrals are performed over a thin shell in (p, lo) space. 
It will be convenient to define the RG so that we inte- 
grate over arbitrary p x , lo but over a shell in p y where 

A > \p y \ > A/6s. The vertex correction becomes 

Su= ~Tf [ D (Py^)SR(Px,Py,^)GL(Px,Py,U}) (35) 

Here Qr/l are the propagators of the right and left mov- 
ing fermions respectively. These propagators include the 
singular frequency dependent self energy discussed in pre- 
vious sections. Note the minus sign in front which comes 
from the different signs with which the gauge field couples 
to the left and right fermions. In contrast at a Pomer- 
anchuk transition the minus sign will be absent. Con- 
sequently what is an enhanced vertex in the gauge field 
problem will become a suppressed vertex at the Pomer- 
anchuk transition. 

We evaluate the integral in Appendix |DJ The result 
takes the form 



Su = 



1 



4tt 2 A 



g(XN, z b ) ln(6) 



where the function g is given by 



g(v, z b ) = I dt ( 

lo \jt+ V t~ +v 2 



Vt Z b 



(36) 



(37) 



where J = t~- For fixed v in the limit that z\, — 2 this 
function has a sensible limit which we demote 32 (v). We 
find 



92{v) 



2['J 2 V 2 + 1) \ 7T7U 



1 ln(jv) (38) 



The vertex correction in Eqn. [551 implies a modified 
scaling exponent for u: 



1 



4tt 2 N 



g{XN, z b ) 



(39) 



As g has a finite limit when zb — > 2, N — > 00 but AA^ 
is finite, this modification is of order l/N. We note that 
it is sufficient to evaluate g at zi, = 2 in this limit. To 
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leading order in e = z\, — 2, jj, we therefore have (using 
A = 27r 2 e appropriate for small e) 



= 1 



1 



4tt 2 A 



g 2 (2ir 2 6N). 



(40) 



Inserting this into Eqns. [521 and [Ml we see that for small 
e, 1/N the power of u> in the 2Kf singularity is 



1 e / 92 {2k\N) 

2 4V TreiV 



(41) 



Since g2(v)/v runs from oo to as v runs from to oo, 
the 2Kf singularity is suppressed over that of the Fermi 
liquid for eN — > oo while it is enhanced when eN — > 
0. The competition in Eqn. [41] between the first and 
second term in parentheses is precisely the competition 
between quasiparticle smearing and Amperean attraction 
described at the beginning of this section. 



The structure of the singularities in the Cooper chan- 
nel will also be modified from that of the Fermi liquid. 
But here the Amperean interaction is repulsive in the 
gauge field case and thus the Cooper singularities will be 
weaker than in the Fermi liquid. On the other hand at 
the nematic quantum critical point the boson mediated 
interaction is attractive in the Cooper channel. Conse- 
quently there will be an enhancement of the pairing ver- 
tex due to boson exchange. This leads to an enhancement 
of the Cooper singularities compared with the Fermi liq- 
uid. A complete discussion of this effect is more compli- 
cated than the 2Kf singularities and will be presented 
elsewhere^ 6 -. Within the large-iV expansion for small 
zi, — 2, in contrast to the 2Kf singularity, there is an 
o(l) correction to the tree level scaling exponent of the 
Cooper vertex which receives contributions from many 
diagrams 36 . Quantitative calculation of this correction 
as well as a study of the interesting interplay between 
superconductivity and critical nematic fluctuations will 
appear elsewhere^ 6 -. 



VI. PERTURB ATIVE FIXED POINT FOR 
FINITE TV, SMALL e AND 2K f SINGULARITIES 

A different controlled limit was previously discussed in 
Ref. Ull. This is obtained by considering finite N, and 
small e where there is a perturbatively accessible renor- 
malization group fermi- liquid fixed point. In this section 
we study the 2Kf singularities at this fixed point. It will 
be useful for our purposes to define a RG scheme that is 
slightly different from Ref. [3lJ. So let us first reproduce 



their main result. Consider the two patch action 

S = S f + S int + S a (42) 
S f = [ d 2 xdT^f sa (d T -isd x -d*) f sa (43) 

sot 



S a 



d 2 xdT^=af sa f s , 



N 



(44) 
(45) 



Compared with the 'minimal' action of Section [TTI above, 
we have set rj — 1 and have reinstated the boson coupling 
e. The ellipses in the last equation refer to other oper- 
ators irrelevant at the RG fixed point to be described 
below (for instance an cu 2 term in the quadratic gauge 
action). When e = 0, the fermions and boson field are 
described by two decoupled Gaussian theories. This is at 
a fixed point under the scaling 

' 7 (46) 



b 

X = b 
, V 

y = Tb 

/V.l/V) = bif(x,y,r) 
a'{x',y',T r ) = b^-^a{x,y,T) 



(47) 
(48) 

(49) 
(50) 



Now turn on a small e ^ 0. By power counting we find 
e' = eb^ (51) 



In differential form (if we let b = 1 + dl) we get 



dl 



(52) 



Thus e 2 is relevant for e > and irrelevant for e < 0. To 
determine the fate of the theory for e > let us study 
the one loop beta function for e 2 . We will define the RG 
by integrating out modes with A > \q y \ > -A^. To order 

e 2 , the fermion self energy is given by the integral 

/ dq y D(q y ,uj')g{p-q ) uj-u}') (53) 

Vb 



E = 



4tt 3 A 



where D is the gauge propagator in the bare action. After 
doing the q x integral the remaining w' integral only gets 
contributions from small u/ so we can replace the gauge 
propagator by its low frequency form ^ ^ b _i . For small 

external frequancy w we get 
£ = - 



ie 2 0J f A dq y 



2tt 2 A 



v 77 



Zh — 1 

Qv 



(54) 



Anticipating that there is a new fixed point when e 2 /N ~ 
o(e), we replace the integrand by its value at z& = 2 to 
get 



ie 2 ui 
'Att 2 N 



In 6 



(55) 
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Thus the inverse fermion propagator takes the form 

„2 



4ir 2 N 



In b) -p x - p 



(56) 



The correction to the fermion-boson vertex (see Fig. [5]) 




FIG. 8: One loop correction to the fermion-boson vertex. 

at order e 2 vanishes as the q x integral has poles only 
on one side of the complex plane. Finally the change 
to the boson propagator is also zero if only modes with 
high \q y \ are integrated out. Thus the only change is in 
lo dependence of the fermion propagator. This may be 
incorporated into a modified scaling 



4tt 2 A 



Info 



■ l- 

LUO 



This implies 



b 1+ < 
x 

b 

y_ 

Vb 



f'(x',y',T r ) = b*+^f(x,y,T) 
a (x ,y' ,t') = fo 5 ~~ + 8^w a(x, y, r) 



(57) 
(58) 

(59) 

(60) 
(61) 

(62) 
(63) 



The modification to the flow of the coupling e is now 
readily obtained to be 

z b -2 „2 

e' = eb 4 s^n (64) 
In differential form this implies the flow equation 
de 2 



ee 
~~2 



dl 2 4tt 2 A 
Thus we indeed find a fixed point when 



el = 2ir 2 Ne 



(65) 



(66) 



Right at the fixed point the scaling equations above are 
identical (to within order e) to those found earlier in Sec- 
tion [HI] and indeed to that expected based on RPA. The 
differences from RPA discussed in earlier sections in the 
fermion propagator come from three loop calculations, 
and hence are not expected to show up till order e 3 . 



The singularities of many physical quantities can be 
usefully calculated within this e expansion and provides 
an alternate controlled limit to the one we have discussed. 
As an illustration let us calculate the boson propagator 
and the fermion self energy. Let the bare value of the 
electric charge at the cut-off scale be eo. The boson prop- 
agator is given by the usual one loop diagram and takes 
the form 



D{q y ,uj) 



1 



e o M 

4tt \q y \ 



W 1 



(67) 



The Landau damping term does not acquire any correc- 
tions from this perturbative answer at least upto the or- 
der to which the RG has been performed. To obtain the 
frequency dependence of the fermion propagator we ex- 
amine the flow of the coefficient of the iuj term calculated 
above. Let us denote this coefficient rj(l) at an RG scale 
I. The calculation above gives the flow equation 



drj 



4ir 2 N 

Combining with the equation for e 2 we obtain 



d (r\e 2 
dl 



2 Ve 



Thus we find 



V (l)e 2 (l) = e\ 



e 2 



(68) 



(69) 



(70) 



where we set r)(l = 0) = 1, and e\ is the bare coupling 
at the cut-off scale. The frequency dependence of the 
fermion self energy is then obtained by setting / = In ^- . 
In the limit oj — » 0, we may set e 2 (/) = e 2 so that 







2tt 2 Nc 



A 2 



Then at order e the self energy becomes 



2ir 2 Ne 



M 1 =sgn(w) 



(71) 



(72) 



Here we have ignored a term A e in the numerator to this 
order in e. 

To compare with the results of previous sections we 
need to set the bare coupling e\ = \. Note in particular 
that the prefactor to the frequency dependence is exactly 
-Ar consistent with the earlier analysis. On the other 
hand to calculate the scaling dimensions of any opera- 
tor directly within this epsilon expansion we need to sit 
right at the fixed point and perturb the theory with that 
operator. The fixed point theory corresponds to setting 
the bare coupling eo = e*. Thus the boson propagator 
and fermion self energy right at the fixed point take the 
forms 



D*(q y ,uj) 



1 



4tt \q v \ + \1y\ 



^sgn(w) 



(73) 
(74) 
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Note that as e — > these fixed point propagators go over 
smoothly into those of the 'decoupled' Gaussian fixed 
point, as indeed they must. 

As it does not seem to be available in the literature, 
let us now calculate the scaling exponent for the 2Kf 
singularity at this o(e) fixed point. As before we add the 
term in Eqn. [25] to the action. By power counting we 
again find 



i = ub 



(75) 



This is modified at leading order of e by the same vertex 
correction diagram as before. Evaluating the integral as 
before we find 



5u 



4ir 2 N 



dq y 



. 77 



dujD*(q y ,uj) 



£*(w) 



% 
(76) 

Here we have written E*(w) = — iT,*(uj). Naively as the 
vertex correction is already order e 2 ~ e we should re- 
place the integrand by its value at e = 0, i.e by the 
fermion and boson propagators at the Gaussian fixed 
point. However at the Gaussian fixed point £ = uj 
and the frequency integral is logarithmically divergent 
at large uj. This signals that the leading order e correc- 
tion to <p u is not analytic in e. To extract it we keep the 
correct boson propagator and fermion self energy at the 
o(e) fixed point calculated in Eqns. [73]- [74] 

Inserting these into the integral for the vertex correc- 
tion, we see that the high-w divergence of the uj integral 
is cutoff by the presence of the Landau damping term in 
the boson propagator. The integral is readily evaluated 
for small e, and we find 



5u = — In | 

2 \weN 



In b. 



(77) 



Thus we get the modified scaling equation 

v! = M 6 1+ i ln (^v). (78) 



so that 



1 + - In ' 



2 \ireN 



(79) 



Inserting into Eqns. I3"2l and IMl we see that the 2Kf sin- 
gularity is enhanced compared to the Fermi liquid at this 
fixed point. Note that this answer for 4> u agrees exactly 
with the result of Section [V] when the limit of e N — >• is 
taken. 



First we notice that at low energies the Fermi surface 
is sharp even though the Landau quasiparticle has been 
destroyed. This is qualitatively the same as in the RPA 
but at least in the nematic case the detailed singular 
structure is modified. Despite this, as argued in many 
previous paper o 3 ' 27 ' 28 this is a compressible state. This 
follows immediately from the general argument in Ap- 
pendix [X] that the compressibility does not receive any 
singular contributions from the low energy scale invari- 
ant fluctuations. Hence this non-fermi liquid state has a 
finite and non-zero compressibility. Actually in the ne- 
matic case changing the chemical potential will also in 
general drive the system away from the critical point. 
This leads to a singular contribution to the ground state 
energy as a function of chemical potential which could 
lead to a singular contribution to the compressibility 23 . 
The safe statement then is that the differential change 
in density in response to a change in location within the 
phase boundary is finite. 

Actually the patch construction implies an even 
stronger result. Consider the susceptibility to a defor- 
mation of the Fermi surface in any angular momentum 
channel (i.e the response to an external field that cou- 
ples to the corresponding shape distortion of the Fermi 
surface). The universal singularities in this quantity are 
obtained by examining this coupling within the patch 
construction. But the patch theory does not know any- 
thing about the angular dependence of the probe field. 
So within each patch this external field couples in the 
same way as an external chemical potential. Conse- 
quently (just as for the contribution to the compress- 
ibility from the low energy density fluctuations) there is 
no singular contribution to the susceptibility in any an- 
gular momentum channel which all stay finite and non- 
zero. For the nematic critical point the only exception is 
the order parameter channel itself (1 — 2 for the d-wave 
nematic). In that case this argument implies that the 
critical behavior of the order parameter susceptibility is 
correctly given by the mean field Hertz answer and re- 
ceives no singular corrections from the fermions. These 
arguments provide a simple explanation of some recent 
results for the nematic critical point obtained through de- 
tailed calculations^. Consider the approach to the quan- 
tum critical point from the symmetric side where there 
is no nematic order. At low energies the corresponding 
metal is described by Fermi liquid theory characterized 
by Landau quasiparticles with an effective mass m* and 
various Landau parameters. On approaching the quan- 
tum critical point standard scaling arguments show that 
the effective mass diverges with an exponent as 



VII. PHYSICAL PICTURE 

In this section we discuss some qualitative aspects 
of the physical picture of the low energy physics, and 
show how we may understand the results of some of 
the detailed calculations presented in previous sections. 



(80) 



In the Fermi liquid phase the compressibility is expressed 
in terms of m* and the Landau parameter as 



K OC 



F9 
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FIG. 9: After any scattering event in the Cooper channel, two 
fermions remain perfectly nested. 



The constancy of k as the critical point is approached im- 
plies that the Landau parameter F® diverges in exactly 
the same way as the effective mass. Applying this rea- 
soning to other angular momentum channels we see that 
the Landau parameters in all angular momentum chan- 
nels (except the order parameter one itself) must diverge 
in the same way as in* so as to give a constant suscepti- 
bility at the critical point. This is exactly the conclusion 
of Ref. QJ. 

In either the gauge model or the nematic critical point 
the patch construction implies that the only singular 
modification from RPA in the charge density response 
happens at the 2Kf wavevectors (modulo the caveat dis- 
cussed above for the q = response in the nematic case). 
As explained in detail in Section [V] apart from the 2Kf 
particle-hole correlations the main modifications from the 
Fermi liquid in the two particle response are in the struc- 
ture of the pairing correlations. Whether the 2Kf and 
pairing correlations are enhanced or not compared to 
the Fermi liquid is largely determined by the Amperean 
rules. In the gauge field case the pair correlations are 
suppressed and the 2Kf potentially enhanced while the 
opposite is true for the nematic critical point. 

Consider now the 1/N 2 calculation of the fermion self 
energy described in Section IIVI The singular contribu- 
tion to the self energy comes the two diagrams shown in 
Fig. [Jj We note that both diagrams may be expressed in 
terms of appropriate two-particle scattering amplitudes. 
Fig.[TT]is a scattering amplitude in the (particle-particle) 
Cooper channel while Fig. [12] is a scattering amplitude 
in the particle-hole 2Kf channel. Based on the physical 
picture dictated by the Amperean rules, we expect that 
in the nematic case the Cooper channel diagram by it- 
self leads to a self energy that is more singular than the 
'bare' terms in the action while the 2Kf diagram by itself 
leads to a singularity that is less singular than the bare 
term. The situation is clearly reversed in the gauge field 
model. This physical picture thus enables us to under- 
stand the signs of the contributions of the two diagrams 
in the calculation. 

Whether the net effect is to produce a singular cor- 
rection to the self energy that dominates over the bare 
one at low momenta and frequencies is determined by 
the competition between the Cooper and 2Kf contribu- 




FIG. 10: After a scattering event in the 2Kf channel, two 
fermions are no longer perfectly nested. 



p p + l p+ q p 




k k — I k — q k 



FIG. 11: The cooper-channel scattering amplitude. 

tions, i.e by the relative magnitude of the contribution 
of the two diagrams. We now argue that the Cooper 
channel always dominates (i.e gives the bigger contribu- 
tion) consistent with the results of the actual calculation. 
To see this consider both two particle scattering ampli- 
tudes when the external lines are right at the Fermi sur- 
face, and initially satisfy the 'nesting' condition. In the 
Cooper channel this means that the total momentum of 
the two incoming particles is zero. In the 2Kf channel 
this means that the incoming particle-hole pair has mo- 
mentum exactly 2Kfi. In the 2Kf channel exchange of 
a boson with momentum q y leads to a new particle-hole 
pair state which no longer satisfies the nesting condi- 
tion (see Fig. [TU)) . Thus after one such scattering event 
the particle-hole pair is less sensitive to the Amperean 
attraction/repulsion mediated by subsequent boson ex- 
change. In contrast in the Cooper channel, exchange of a 
boson with momentum q y preserves the nesting condition 
for the resulting particle-particle pair (see Fig. . Thus 
they are able to continue to reap the benefits of the Am- 
perean interaction in subsequent scattering events. This 
explains why the Cooper channel always dominates over 
the 2Kf channel. This difference between the kinemat- 



p P + q P + l P 




k k + q k + I k 



FIG. 12: The particle-hole channel scattering amplitude. 
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ics of the Cooper and 2Kf scattering channels is clearly 
present only if the Fermi surface is curved. Thus we 
expect that in the artificial limit where we ignore the 
curvature term in the fermion Greens function the two 
diagrams will have the same magnitude and hence will 
cancel. Examining the relevant integrals shows that this 
is exactly what happens. 

Let us now discuss one important physical consequence 
of these results. At the nematic critical point, the sin- 
gular structure of the fermion Greens function is modi- 
fied from RPA at order 1/N 2 . Specifically the fermion 
Greens function satisfies the scaling form of Eqn. [7] with 
the fermionic dynamical exponent z = % and the expo- 
nent a = 1 — f]f with r\j given in Eqn. [21] Note that r]f 
is positive. This represents a modification of RPA which 
has r/f — 0. As pointed out in Ref. |32| . this shows up very 
directly in the tunneling density of states N(uj) defined 
through 

*(«>) = J 0y A ^^) ( 82 ) 

where the single particle spectra function A(K, u) is de- 
fined in the usual manner: 

A(K,u) = lmG(K,iw^uj + iO + ) (83) 

Here K is the full momentum (not linearized near the 
Fermi surface). Apart from being a potentially direct 
measure of the deviation from RPA, study of N(u) also 
provides insight into the sign of rjf and some rational- 
ization for why it is non-zero in the first place. Singular 
contributions to N(uj) come from momenta in the vicin- 
ity of the Fermi surface. The two dimensional momentum 
integral may then be separated into an angular integral 
over the Fermi surface and a radial integral over just 
the component of the momentum normal to the Fermi 
surface. The former just contributes an overall constant 
prefactor. The latter may be directly evaluated to obtain 
the result advertised in the introduction 

N(u)~\u\¥ (84) 

Thus at the nematic critical point there is a power law 
suppression of the local single particle density of states. 
This suppression is of course rather natural if we remem- 
ber that superconducting fluctuations are enhanced at 
the nematic critical point. Thus the sign of r\f may be 
qualitatively understood. Further the enhanced super- 
conducting fluctuations make it plausible that there be 
some effect on the density of states unlike what happens 
in the RPA. 



VIII. TOWARDS A PHASE DIAGRAM 

We now turn to the question of what happens for gen- 
eral Zb,N. The preceding sections show that so long as 



Zb — 2 is small, the theory can be controlled for any N. 
What happens if Zb — 2 is not small? As discussed in the 
introduction for Zb = 3 the possibility of using large- A as 
a control parameter has been studied in detail recently 
and several difficulties have been pointed out. Here we 
suggest an interpretation of these difficulties. Consider 
first the nematic critical point. The action for the order 
parameter fluctuations within the t wo p atch theory was 
calculated to three loop order in Ref. [23|. For fluctuations 
at zero frequency they found 

J^N (l - cy/ltj \k y \ 2 \a(k,u = 0)| 2 (85) 

The second term comes from the three loop polarizabil- 
ity of the fermions. The appearance of the extra factor 
of y/N in the loop calculation raises concerns over the 
existence of a sensible large- N limit. It is currently not 
known what the structure of higher loop terms are - for 
instance whether same or even higher powers of N are 
generated by higher loop contributions. If we take the 
three loop answer at face value then as c > 0, for large 
enough N, the coefficient of k 2 , becomes negative. This 
signals an instability towards ordering at non-zero mo- 
mentum. In particular this means that the original as- 
sumption of a direct second order nematic transition is 
not correct, and the nematic transition will be preempted 
by the appearance of density wave order. 

Can this conclusion be changed by higher loop dia- 
grams? One possibility is that higher loop diagrams 
change the sign of the coefficient c. But then a differ- 
ent instability will likely set in. For instance consider 
diagrams with the same structure as those in Fig. [5] but 
with arbitrary number of boson lines connecting the right 
and left moving fermions. Whatever the sign of the sum 
of diagrams of this sort, so long as it has a higher power 
of N than the one-loop one, there will be an instabil- 
ity. If the sign is negative (as in the three loop calcula- 
tion), then there is an instability where the boson likes 
to order at non-zero wave vector. If the sign is positive 
then we consider the response to an external probe field 
that couples with opposite sign to the two patches. A 
concrete example at the nematic critical point is just an 
external electromagnetic gauge field, i.e we consider the 
'diamagnetic' response to a static external magnetic field. 
This changes the sign of the external vertices in diagrams 
like Fig. [5] without changing the sign of the internal ver- 
tices. The 'diamagnetic' response to a static magnetic 
field then has the opposite sign from that of an ordinary 
metal. This signals an instability toward spontaneous 
flux formation, i.e the system will likely develop a state 
associated with spontaneous circulating currents. 

The only remaining possibility is that higher loop dia- 
grams exactly cancel the offending y/N term found in the 
three loop calculation. While we cannot rule this out we 
can provide a suggestive argument against this possibil- 
ity by examining the limit of small Zb — 2. In this limit, 
the static boson polarizability may formally be written 
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as a series 

n( ky , u = o) = |jy«*-i + £ /n( ffi zt) ) (86) 

The leading n — 1 term was calculated in Appendix [C] 
Let us assume that the functions f n {x; Zb) all have finite 
limits when Zb — > 2: 

lim f n (x,Zb) = F n (x) (87) 

This is explicitly seen to be true for n = 1, and we assume 
it holds for arbitrary n. Then successive terms in the 
series above are down by powers of 1 /N for finite non-zero 
XN. So for large- N there is no instability. The instability 
potentially happens when XN becomes large enough that 

2 

the n = 1 term is comparable to 1, i.e when XN ~ . 
If any higher order term, say the nth one, is to have the 
same power of N in the large AiV limit (while still keeping 
e small), then F n (x) ~ x n+ ^ for large x. But then its 
coefficient has a high power of A ~ e. Therefore any 
cancelation of the dangerous three loop term by higher 
loop diagrams cannot in general happen for arbitrary Zb- 
This makes it rather likely that there is an instability for 
arbitrary Zb not too small. 

For Zb approaching 2 comparison of the three loop term 
with the leading one loop term suggests that the insta- 
bility happens when e« ~ -j^. This leads to a 'phase 
boundary' between the unstable and stable regions that 
comes in with infinite slope in the e, plane. Thus we 
propose the phase diagram shown in Fig. Q] or Fig. [3] 
Through out the unstable region there is no direct ne- 
matic transition, and it is always preempted by instabil- 
ity toward a different broken symmetry. It is not clear 
whether the unstable region encompasses the all impor- 
tant point N — 2, Zb — 3, i.e whether Fig. Q] or Fig. [2] 
is realized. If however this point belongs to the stable 
region, i.e Fig. [1] applies and there is a direct second or- 
der nematic transition, then we have no choice but to 
access it from the small e region (either by combining 
with large- N, or by the perturbative RG for small N). 
If Fig. [2] is realized on the other hand we may still hope 
that our expansion captures the physics at temperatures 
above the instability. 

Similar considerations apply to the gauge field model. 
There at Zb = 3, the gauge polarizability acquires only 
a positive o(y/~N) correction at three loop level. How- 
ever at the same order if we consider the response to an 
external probe that couples with the same sign to both 
right and left movers, then the sign of the three loop 
response is reversed. In arguing for an instability, it is 
important that the bare zero-momentum susceptibility 
vanish, since the result Eqn. [35] applies in the scaling 
regime, which requires \k y \ <C N 3/2 V if the bare suscep- 
tibility were nonzero and TV- independent, the putative 



instability would occur at k y outside this regime 2 . A 
suitable choice may be a perturbation of both the volume 
of the system and the chemical potential, preserving the 
particle density. If the three loop calculation were the 
full story this would again signal an instability toward a 
state which spontaneously orders at non-zero momentum 
and hence breaks translation symmetry. More generally, 
when higher loop terms are included the situation is sim- 
ilar to the discussion above for the nematic model. Con- 
sequently we suggest that the gauge field model is also 
unstable toward a state with some broken symmetry at 
sufficiently large- N when Zb is sufficiently different from 
2. Thus once again this proposal would imply that if the 
gauge model at N = 2, Zb = 3 is stable then we have no 
choice but to access it as we have done from the small e 
region. 



IX. DISCUSSION 

In this concluding section we consider the implications 
of our results for some specific systems, and for the gen- 
eral theory of non-fermi liquid metals. 

An important and topical realization of the gauge 
model is to the theory of gapless quantum spin liquids 
where a gapless Fermi surface of charge neutral spin- 
1/2 fermionic spinons is coupled to a gapless fluctuat- 
ing U(l) gauge field. Note that in this example the spin 
liquid is a non-fermi liquid metal for spin transport but 
is an insulator for electrical transport. Such a state has 
been proposed^* 2 - to describe the intermediate tempera- 
ture scale physics of the layered organic Mott insulators 
k - (ET) 2 Cu 2 (CN) 3 and EtMe 3 Sb[Pd(dmit) 2 ] 2 . Our 
controlled calculations merely confirm the correctness of 
several key results from RPA that are directly relevant 
to experiments - for instance the scaling structure of the 
low energy theory implies that the specific heat follows 
the familiar RPA result C v ~ at low temperature 
T. The more important contribution of the present pa- 
per to the theory of such a spin liquid is the controlled 
calculation of the structure of the 2Kf spin correlations. 
Detecting these in experiments would be an interesting 
way to 'measure' the spinon Fermi surface (see Ref. |37J 
for a proposal) . Our results also set the stage for an anal- 
ysis of phase transitions from the spinon Fermi surface 
state to various proximate phases with spinon pairing or 
other 'order' that may be relevant to describing the very 
low temperature physics of the organics. 

The model of a fermi surface coupled to a gauge field 
also describes algebraic charge liquid metals^&i and the 
related d-wave Bose metals^. An essential difference with 
the particular gauge model studied in this paper is that 
there are two species of fermions that couple with op- 



2 We thank Max Metlitski for helpful correspondence on this 
point. 
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posite gauge charges to the same fluctuating U(l) gauge 
field. The Amperean rules are therefore different and this 
will lead to some differences in the results. These can be 
straightforwardly handled within our expansion. Simi- 
larly our methods are readily generalized to provide con- 
trolled expansions for various Pomeranchuk transitions 
other than the nematic example considered in detail in 
this paper. 

It is interesting to contrast the quantum critical metal 
at these Pomeranchuk transitions with other examples 
of non-fermi liquid metals. One other set of examples is 
provided by Mott-like quantum phase transitions where 
an entire Fermi surface disappears continuously. Apart 
from continuous Mott metal-insulator transitions, these 
are thought to describe non-fermi liquid physics in heavy 
fermion metals near the onset of magnetic long range or- 
der. Ref. |32| argued that such Mott-like transitions will 
be characterized b the presence of a sharp critical Fermi 
surface but without a sharp Landau quasiparticle. A crit- 
ical fermi surface is also a feature of a Pomeranchuk tran- 
sition if it is second order. However we might expect that 
the destruction of the Landau quasiparticle is more severe 
at the Mott-like transitions. Indeed the explicit calcula- 
tion in Ref. HO for a continuous Mott transition found the 
exponent value a = — r\ (where 77 is the anomalous expo- 
nent of the boson field at the 3D XY fixed point, and is 
known to be small and positive). This corresponds to a 
large fermion anomalous dimension rjf = 1— a = 1+?/. In 
contrast at the nematic critical point, the fermion anoma- 
lous dimension is small. Within RPA it is simply while 
the three loop calculation of Ref. [23| as well as the con- 
trolled estimate presented in this paper give non-zero but 
small values. The largeness of rjf is a partial measure of 
the extent to which the quasiparticle is smeared (as is 
exemplified by the suppression of the tunneling density 
of states). 

It is instructive to compare the non- Fermi liquids stud- 
ied in this paper with those discovered recently using 
holographic duality ^ 41 ' 43 , which may be understood 
heuristically as arising from a Fermi surface coupled with 
some bath of critical fluctuations with an infinite dynam- 
ical exponent^ - — . It would be interesting to find explicit 
field theoretical models with this feature. In the models 
studied here, large Zb > 3 is unstable because a \q\ 2 term 
in the boson inverse propagator will always be generated 
by short distance fluctuations and will eventually domi- 
nate. 

In summary in this paper we have developed a con- 
trolled and systematic approach to calculating the uni- 
versal properties of a non-fermi liquid metal that arises 
when a gapless Fermi surface is coupled to a fluctuat- 
ing gapless boson field. We illustrated our approach by 
studying spinon fermi surface spin liquids, and quantum 
critical metals near an electronic nematic transition in 
some detail. Our approach readily lends itself to the 
study of various closely related problems. We leave the 
exploration of these to the future. 
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Appendix A: Universality and the patch 
construction 

In this paper we have discussed a scaling theory which 
focuses on the interactions of fermions near the Fermi 
surface via bosons of small frequency and small momen- 
tum. The kinematics of these bosons allows us to re- 
strict attention to one patch of the Fermi surface and 
its antipode. In particular the universal singularities in 
the low energy physics is correctly captured by break- 
ing up the full Fermi surface into patches (and their an- 
tipodes), and studying the theory patch by patch. In 
this Appendix we briefly discuss some subtle points as- 
sociated with the patch construction and the subsequent 
treatment of the theory for any given pair of antipodal 
patches. As a bonus we will show that the compress- 
ibility has no singular contributions coming from the low 
energy scale invariant fluctuations. 

Consider the following microscopic Lagrangian for the 
gauge field problem 



C 



uv 



1 

2m 



-*V + -La) -// ) / . (Al) 



This Lagrangian is gauge-invariant, and this forbids a 
mass term for the gauge field a in the effective action. 

Next consider the low-energy description which focuses 
on the modes near two antipodal patches of Fermi sur- 
face (with normal x, without loss of generality), / ~ 
f R e ikpX + f L e- ikFX : 



-'patch 



ft[d T - ivpd x 



d 2 

2m 



VF 



a fa 



+ fUd T + iv F d x - - ^La ) f L (A2) 



2m 



N 



and vf = kp/m. Here a is the ^-component of the gauge 
field which couples strongly to this pair of patches. For 
other patches with normal n(9) at an angle 9 to the x- 
axis, it is the component a-h{9) that will couple strongly. 
The full theory is obtained by summing over all patches 
and adding together the 'diamagnetic' term 



C 



dia 



2mA 



(A3) 
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By itself the patch action appears to respect an 'emer- 
gent' gauge symmetry which acts by 

a^a + dX, fn-te'jkfit, h^e^h (A4) 

with A = X(x). However, we observe that the gauge field 
a couples in this two-patch theory to (vp times) the axial 
current /Jj/r — /J-/l, which is "anomalous", as we now 
review. To diagnose this anomaly, consider the coefficient 
of a 2 in the effective action resulting from integrating out 
the fcrmions. In the patch approximation, the numbers of 
left- and right-moving fermions are separately conserved. 
Turning on the gauge field a violates this conservation 
since the gauge field couples like a chemical potential 
with opposite sign on the two sides. The change in the 
'chiral density' of fermions is then the density of states at 
the Fermi surface times the effective chemical potential 
change, . 

However in the patch theory, the density of states at 
the Fermi surface is ill-defined. It is apparent from its 
high-energy origin as a theory with a finite Fermi surface 
that this description is only valid up to some maximum 
deviation of the momentum from the middle of the patch. 
Let the patch size (which equals the cut-off for q y ) be 
denoted A y = KpA8 where A6 is the angular extent of 
the patch. The density of states in each patch is then 
N ^i^ e = E j£ 1 - This apparently implies that the one- 
loop gauge field polarizability 

n(£,w = o) = (o(M)o(-M))i_ Wp (A5) 

depicted in Fig. Stakes a non-zero value as k — > 0: 

n(fc^0,w = 0) = (A6) 

Z7T Z 

Note that this is the contribution from both patches to 

2 

the fermion polarizability to the coefficient of — ^a 2 in 
the one loop euclidean effective action for the gauge field. 
Naively this seems to be a problem as it violates the fake 
gauge invariance (|A4[) ; more problematically, it also vio- 
lates the real microscopic gauge invariance. The resolu- 
tion is that the microscopic gauge invariance is obtained 
only when the diamagnetic term is also included in the 
effective action. Indeed if we sum over all patches the 
contribution from the polarizability is 

- / d9— F -(a-nf = -—^a 2 (A7) 
J 47r 47rm 

As the density satisfies po = N ^ F ■ we see that the non- 
zero fermion polarizability at k —¥ 0, uj — exactly can- 
cels the diamagnetic term as required by gauge invari- 
ance. 

The answer above for the fermion polarizability may 
be reproduced formally by considering the integral 

n(fc,w = o)cx f g(k+ P )g(p) (as) 

Jp 



Actually this integral is ill-defined at short distances and 
depends on the order of integration. The origins of the 
patch theory from the original full Fermi surface means 
that we need to impose a cut-off on the x and y momenta 
but not necessarily on the frequency. Doing the frequency 
integral by contour integration we reproduce (in the limit 
k — > 0) the non-zero constant obtained above through a 
physical argument. Note that imposing a hard cut-off on 
p x violates the fake gauge invariance of Eqn. IA4I Thus 
a careful formulation of the patch theory that is faithful 
to its microscopic origins requires imposing a cut-off on 
the fermion momenta, and taking the frequency cut-off 
to infinity first before sending the momenta cut-offs to 
infinity to define the scaling limit. 

However as in any scaling theory we expect that uni- 
versal low energy singularities are actually insensitive to 
how the theory is regularized at short distances (i.e in- 
dependent of the "short distance completion" of the the- 
ory). Thus for the purpose of calculating the universal 
singularities we are free to choose any regularization of 
the patch theory. A convenient choice (and the one used 
in this paper) is to define the scaling theory by send- 
ing the momenta cut-offs to infinity first, and then the 
frequency cut-offs. Then the polarization integral may 
be done by first integrating over p x . The result then 
vanishes as both poles of p x lie on the same side in the 
complex plane. Thus in this regularization of the patch 
theory there is no need to include a diamagnetic term to 
maintain gauge invariance. Indeed with this regulariza- 
tion the fake gauge invariance is no longer fake and is a 
real property of the universal scaling theory. 

Though this choice for defining the scaling theory will 
reproduce all the universal singularities there is no guar- 
antee that it will correctly reproduce non-universal ones 
as it is not faithful to the original microscopic action 
with the full Fermi surface. A good example to illus- 
trate this is the fermion compressibility, i. e the response 
in the density of the system to a change of chemical po- 
tential. Within the patch construction the chemical po- 
tential couples to f R fn + f\fh- This is identical to the 
coupling of the nematic order parameter. In the single 
patch theory, the chemical potential term is thus identical 
to the gauge coupling term. The contribution to the com- 
pressibility is thus given by the same n(fc — > 0,o> = 0) 
discussed above. With the regularization actually em- 
ployed in this paper (which dictates us to do the p x inte- 
gral first) this is zero as is indeed required by the gauge 
invariance of Eqn. IA4I In the other regularization (dic- 
tated by the original microscopic situation) we do the 
frequency integral first and get a non-zero answer. 

What is the interpretation of the zero answer that do- 
ing the p x integral first produces? Clearly it means that 
there is no singular contribution to the compressibility 
in the low energy theory. The non-zero answer found 
in the other regularization must be viewed as a smooth 
non-singular background that is not correctly described 
by the scale invariant patch theory. The validity of Eqns. 
IA4I in the scaling limit employed in the paper thus en- 
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sures that singular contributions to the compressibility 
vanish to all orders of perturbation theory in the one 
patch theory. Actually the same result is also true in the 
two patch theory. To see this note that with the short 
distance completion we have chosen the two patch action 
in the presence of an external potential that couples to 
JrIr + /i/i enjoys a low energy gauge invariance sim- 
ilar to Eqn. IA4I but where we rotate the phases of right 
and left movers with opposite phases ±7^- This gauge 
invariance ensures that the contribution to the fermion 
compressibility vanishes even in the two patch theory. 
For the original microscopic model this then implies that 
there is no singular contribution to the background non- 
zero compressibility. 

In Appendix [Ej we explicitly check that short-distance 
sensitivity (similar to what happened in the polarizability 
integral) does not arise in some of the universal quantities 
of interest, and the answers are insensitive to the order 
of integration. 



Appendix B: Example diagram - single patch 

Consider the 3-loop self-energy contribution shown in 
Fig. [T2] with the external Fermion on the Fermi-surface, 



we hnd 



("+" and "— " for the right and 



left patch, respectively) for Zb — > 2. 




k + q 



FIG. 13: Three-loop fermion self-energy containing fermions 
from a single patch. 
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' v qlydbjqdu>isga{q v — l y ) 
D(q)D(l). (B4) 



One still needs to verify that the remaining integral is fi- 
nite. In Ref I22I it has been shown that all planar diagrams 
in the single-patch theory are UV-finite, so it remains to 
verify IR-fmitcness. For I finite and q — >• or q — >• I 
this is clear. For q, I — > the integrand diverges as q~ 2 
which is cancelled by the remaining integrations. Thus 
the 3-loop self energy is indeed of order N~ 2 as long as 
XN is of order unity. 



Appendix C: Self energies at three loops 

Ref. [H identified the important diagrams that con- 
tribute to the boson and fermion self energies at three 
loop level, and evaluated them in appropriate limits when 
2h = 3. Here we briefly sketch the modifications for gen- 
eral Zb close to 2. Consider the diagrams in Fig. [5] that 
determine the boson self energy at three loops. For Zb f=s 2 
we find 
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In the physical limit 7A = 
integral numerically to find 



s) 2 + s 2 
(CI) 
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n = o.i06|fc y | 2:fa - 1 . 



2 we evaluate the 



(C2) 



(Bl) Now consider the diagrams for the three loop fermion 
self-energy contributions depicted in Fig. [7] 

xGR/L(k + q)g R/L (k + i)g R/L (p + q)g R/L (p + i), They are given by 



where 
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Begin by integrating out p — (ui p ,p) 
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and dl = 



[_ J_ J D(q)D{p)D{q-p) (C3) 

* J UJi,l J LO q ,q J Wp,P 

x g R (k ± P )g R (k ± q )g L {i)g L (i + q)g L {i + p). 

The integrals may be evaluated following Ref. [23|, the 
results are, to leading order in 
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where / is a non-singular function independent of Zb and 
N. The l x and q x integrations are straightforward and 



where 
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The function J(XN) which appears in Eq. (|24| is the 
sum of the singular contributions at Zf, = 2, i.e. 



J(A.Y) (^] (/',(()) -r,(A.V)) (C7i 



(note that the prefactor is unity at z\, = 3, N = 2). A 
numerical estimate of this function is shown in Fig. 1141 
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zero. The computation is straightforward: 
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The corresponding expression at the perturbative (small 
e) fixed point is given by 
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FIG. 14: Numerical integral determining the fermion anoma- 
lous dimension r\; . 



The integral remains finite if we take the limit e — > 
in the first term in the integrand, so we can evaluate it 
analytically. We obtain to leading order in e 



Appendix D: Vertex integral 
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FIG. 15: One loop correction to the 2fcF-vertex. The dotted 
line denotes a 2/cf density fluctuation. 

We want to evaluate the integrals shown in Fig. IT5| 
where all external momenta and fequencies are put to 



Appendix E: Order of integration 

Consider as a concrete example the three-loop self en- 
ergy depicted in Fig. [7] We want to keep a cut-off on 
they y-components of momenta and verify that we can 
interchange the order in which we perform the integra- 
tion over frequencies and momentum cc-components. As 
the only possible divergencies arise from the UV, we can 
set all y-components as well as external momenta and 
frequencies to zero (the y components in the numerators 
of the gauge propagators are set to unity). Recall that 
any two integrals are independent of the order of integra- 
tion, if the double integral is absolutely convergent. So 
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consider 

/ iDugDug-uiiD^GpGqGlGp+qGp+ll (El) 

Since each of these integrals is individually convergent, 
the original integrations can be performed in any order. 



This analysis applies to most of the diagrams shown here 
and is a result of the Landau-damping form of the gauge 
propagator, thus it does not extend to diagrams calcu- 
lated at the small e pertubative fixed point. In our pro- 
posed expansion N — > oo, XN = o(l), of the diagrams 
shown in this paper only the one-loop boson self-energy 
(Fig. [3]) is sensitive to the order of integration. 
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